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1 Introduction 

The spinless Salpeter equation is a simple relativistic version of the Schrodinger 
equation which can be obtained, with some approximations, from the covariant 
Bethe-Salpeter equation [1,2] and takes the form {h = c= 1) 



where m is the mass of the particle and M is the mass of the eigenstate (M = 
m + E, E is the binding energy). We restrict our attention to time component 
vector potentials. This equation is generally used when kinetic relativistic 
effects cannot be neglected and when the particles under consideration are 
bosons or when the spin of the particles is neglected or is only taken into 
account via spin-dependent interactions. Despite its apparent complexity, this 
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equation is often preferred to the Klein-Gordon equation. The equation (1) 
appears, for example, in mesons and baryons spectroscopy in the context of 
potential models (see for example [3,4,5,6,7,8,9,10], see also [11]). 

Due to the pscudo-diffcrcntial nature of the kinetic energy operator, few ex- 
act results arc known about this equation. Most of these results have been 
obtained for a Coulomb potential (for example, upper and lower bounds on 
energy levels) [12,13,14,15,16]. Recently, upper and lower limits on energy lev- 
els have been obtained for some other particular interactions [17,18]. 

Conversely to the Schrodinger equation, for which a fairly large number of re- 
sults giving both upper and lower limits on the number of bound states can be 
found in the literature (see for example [19,20,21,22,23,24,25,26,27,28,29,30,31,32,33,34,35]), 
only one result, concerning the total number of bound states, is known for the 
spinlcss Salpeter equation [36]. After recalling in Section 2 a general method 
to obtain upper limits on the number of bound states due to Birman [20] and 
Schwinger [21], we derive such limits for the spinless Salpeter equation in Sec- 
tion 3. The Section 4 is devoted to test these hmits. In Section 5, we show how 
a simple modification of the limits found in Section 3 leads to upper limits on 
the number of bound states that lie below a given energy. 



2 The Birman-Schwinger method 

Birman [20] and Schwinger [21] have shown how to obtain an upper limit on 
the number of bound states once the Green function of the kinetic energy 
operator of a wave equation is known. In this section, we recall briefly the 
main hue of the method; for more details see the original articles [20,21]. 

Let T (p^) be a general kinetic energy operator and let 




(2) 



be the wave equation, in three dimensions, that satisfy the wave function ^(r) 
(eigenstatcs) and where E is the energy (eigenvalues). Let G{A) be the Green 
function of T (p^): 




where A = r — r', A = | A| and (5'^(x) is the Dirac function. We can write (2), 
with E — 0, using this Green function as 





Since the purpose of the method is to obtain an upper limit on the num- 
ber of bound states, we can replace V{r) by — |y~(r)| where V~{r) is the 
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negative part of the potential obtained by setting the positive part of the po- 
tential eqnal to zero. Indeed, a decrease of the potential in some region must 
lower the energies of the bound states and therefore cannot lessen their num- 
ber. Moreover, we introduce the parameter < A < 1 by the substitution 
|y~(r)| — > A|y~(r)|. As A increases from 0, we reach a critical value, Ai, at 
which a bound state first appears with a vanishing binding energy, E = 0. 
With further growth of A, the energy of this state decreases until we reach a 
second critical value, A2, at which a second bound state appears and so on. 
When A has attained the value unity and, Aat < 1 < Xn+i, there are N bound 
states. 

We now introduce, to obtain a symmetrical kernel, a new wave function as 

$(r) = |F-(r)|V2^(r). (5) 

The equation (4) becomes 

<l>(r) = \Jdr' K{r, r') <l>(r'), (6) 

where ^(r, r') is given by 

K{r,r') = \V-{r)\^/^G{A) \V-{r')\^/\ (7) 

If the kernel is positive, we have 0<Ai<A2<---<AAr<l and < A^ < 00 
{\k denotes each eigenvalue of (6)). It is well known that the trace of the 
iterated kernels equals the sum of the eigenvalues of the integral equation (6) 
as follow 

00 -1 „ 

T.7V^= drK(-\r,r), (8) 



where the iterated kernel X^"^ (s, t) is given by 

X (") (s, t) = Jdu K{s, u) K'^''-^^ (u, t) , (9a) 

with 

xW(s,t) = K(s,t), (9b) 
and n — 1, 2, . . .. Now it is plain that the following inequalities hold 

CO 1 A'' 1 

E TtV > E > (10) 

fri (A.)" - fri (Afe)'^ 

where is the number of bound states. Prom (8), (9) and (10) we find that 
an upper limit on the total number of bound states of the wave equation (2) 
is given by 

N < J cirA:(")(r,r). (11) 
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3 Application to the spinless Salpeter equation 



3. 1 Upper limit on the total number of hound states 



To obtain an upper limit on the total number of bound states of the spinless 
Salpeter equation, we need to calculate the Green function of the kinetic energy 
operator. Similar calculations have already been performed previously [37,38]. 
In contrast to calculations found in [38] , we need here to calculate the Green 
function of the following operator 



T (p2) = yp^Tm2 - m. (12) 
This is done by performing the integral 

GKA)^-L,/dp ;-P<-'P-^) . (13) 
(ztt)'^ J VP + — m 



We find that 

' m 



with 



2 

1 + -F{mA) 

TT 



F{y) = I T^i(^) + 2' (1^^) 

where K^{y) is a modified Bessel function (see for example [39, p. 374]). The 
zero-energy spinless Salpeter equation takes thus the form of the following 
integral equation 

*(r) = J dr'G{m,A)\V-{r')\^{r'), (15) 

with G{m, A) given by (14), and where we have replaced the potential by its 
negative part. This integral equation (15) can be written with a symmetrical 
kernel: 

$(r) = J dr' K{m,r,r')<^{r'), (16a) 

where 

K{m,r,r') = \V-{r)\^/^G{m, A)\V-{r')\^/\ (16b) 

and $(r) given by (5). We can now use the Birman-Schwinger method, and 
in particular the relation (11), with the integral equation (16) to obtain the 
following upper limit on the total number of bound states of the spinless 
Salpeter equation 

N dv^...dvr, |y-(ri)| . . . |y-(r„)| ^(m, A12) . ..G{m, A„i), (17) 



47rA 



if(mA), (14a) 
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where Ajj = |rj — rj| and n > 4 (the integral diverges for smaller values 
of n). We have introduced in (17) the parameter a which takes the value 1 
respectively 2 for one respectively two (identical) particle problems. 

Now, we need to calculate the Green function (14) in a closed form, that is to 
say, to compute F{y), see (14b). An integration by part leads to 

Fiy) = K^{y) + ^- H dzK^^z). (18) 

Z Jy 

To obtain an upper limit for the number of bound states, a majorization of 
the kernel (16b) is enough. Since the Bessel function Kq{z) is positive for 
< ^ < cx), and that its integration between and oo is equal to 7r/2 (see 
[39, p. 486]), the integral in (18) is not only positive but also small compared 
to the other terms of this equation. Indeed, in the region [y ^ 0) where this 
integral takes its maximal value (7r/2), this quantity is still small compared to 
the value taken by the singular Bessel function Ki{y). Thus the majorization 
obtained by replacing the integral in (18) by 0, namely 

F{y)<K^{y) + n/2, (19) 

should not spoil too much the upper limit as one can verify with an examina- 
tion of the results of the tests performed in Section 4. Another majorization, 
which leads to a simpler kernel, is obtained by replacing the singular Bessel 
function Ki{y) by 1/y. This additional approximation is however exact in the 
case of a vanishing mass (m = 0) since mKi{my) — 1/y in this limit. The 
kernel is then given by 



The upper limit (17) can then be used with either G^-^^m, A^) or G^'^\m, Aij). 
Of course, using the function G*^^-* {x, y) yields more stringent results than those 
obtained by using the function G^'^\x,y). 

We will discuss more specifically the case of a central potential in the Section 
3.2 but we already write here the upper limit (17), using the kernel G^'^\x^ y), 
for this class of potentials. In this way, we introduce some quantities which 
will be useful later. Integrating over angular variables, the limit (17) reads 



Ar< ^(2^/+l) / dn...drn\V-{n)\...\V-{rn) 



X 



xA^(m,ri,r2) . . . A^(m, r„, ri), (21a) 
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with 



and 



A^{m,x,y) = --^^xya^,{m,x,y), 
a Zu + i 



(21b) 



m ,^-(.+1)^ 1 r(i/ + l) {xyY 



F I/ + 1,1/ + 1,2(1/ + 1), 



27r3/2 r(i/ + 1/2) {x + |/)2(^+l) 
4x?/ \ 



{x + yYl ' 



(21c) 



where r< = min[a;, y], r> = max[a;, j/] and F{a,b,c,z) is the hypergeometric 
function. 

A simpler Tipper limit can be obtained with additional approximations in the 
ultrarelativistic case {m — 0). Using n times {n > 4) the Holder inequality we 
obtain 



roo poo 

N<B(n,p,p') / drr'(f-')/f |y-(r)| / dr\V'(r)\PP 
Uo J Uo 

' [/ 

Jo 



POO 

Jo 



drr"^-^ \V-ir) 



(n-3)/p 



i/=0 



(22a) 



with the constant B{n,p,p') given by 

oo 

B{n,p,p') = E(2^+ 1) 1))]""'^^(^^,Pp7(p' - 1)), (22b) 

and 



C{u, q) 



1 r(;/+l) 
av^22^+2r(t/ + 3/2) 

jj^' dx (1 + l/x^) (z"+iF(i/ + 1, 1/ + 1, 2(1/ + 1), z))' 



, (22c) 



with 2; = 4a;/(l+a;)^,p> 1 and p' > 1. All the complexity of the problem is 
now located in the calculation of the constant B{n,p,p'). Analytical calcula- 
tions for small values of u and numerical investigations up to z/ = 100 seems 
to prove that for all values of v we have 



C(i/,2) = l/{aV2i7TT), 



(23a) 



C(i/,3)<l/(a(2z/ + l)i/3). 



(23b) 
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In this case, we find that B{n. 2, 3) < (1 — 2~'^)([w). with w = (3n — 7)/6 and 
(^{x) is the Riemann Zeta function. This implies n > 5 to obtain nontrivial 
results and in particular B{5, 2, 3) < 2.172. 

Daubechies has obtained the following simple upper limit on the total number 
of bound states [36] 

N <K J dr [\V'{r)\{\V-{r)\ + 2m)f\ (24) 

with K = 0.239 for arbitrary values of m and K = 0.103 for ni = 0. This 
inequality shows that N grows with strength of the potential, g, at most as 
g^. The upper limits obtained in this section behave as g"", with n > 4, and 
should not yield very cogent results for large value of g when the potential 
possesses many bound states. Nevertheless, in the Section 3.3, we obtain an 
upper limit on the total number of bound states for central potentials, and 
m — 0, which features the correct dependency on g. 

3.2 Upper limit on the number of i-wave bound states 

To obtain an upper limit on the number of i-wecve bound states of the spinless 
Salpeter equation, we need to derive the radial version of the integral equa- 
tion (16). This can be achieved by performing similar calculations than those 
done in Ref. [38, p. 2255-2257]. Integration over angular variables reduces the 
integral equation (15) to the following one-dimensional integral equation 

ue{r)= dr'Ge{m,r,r')\V-{r')\ue{r'), (25) 
Jo 

with 

GAm,r,r') = ^ P dO' sinO' ^^^^ PeicosO'), (26) 
2 Jo A 

where 'U^(r) is the radial wave function, \['(r) = {u£{r) /r)Y£m{r) and H{x) is 

defined by (14a). As we pointed out in Section 2, the integral equation (25) 

can be written with a symmetrical kernel provided we introduce a new wave 

function 

Mr)^\V-{r)\'/'ue{r). (27) 
This change of function leads to the following integral equation 

roo 

Mr)= dr'K^{m,r,r')Mr'), (28) 
Jo 

with 

Ki{m,r,r') = \V- {r)\^/^ Ge{m,r,r') \V-{r')\^/^. (29) 

We can now again use the Birman-Schwinger method, see Section 2, with the 
integral equation (28) to obtain an upper limit on the number of i-wave bound 
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states. The relation (11) together with the majorization (19) lead to 

POO 

Ne< / dri. . .drn\V~{ri)\ . . .\V~{rn)\Ti{m,ri,r2) . . .Te{m,rn,ri] 
Jo 



with 



and 



aTi{m,r,r') = -Qi{'m,r,r') + Si{m,r,r'), 

TT 
rr+r' 

Qeim, r, r') dy Ki{my) Pe 

J\r—r'\ 



^2 I y»/2 y"^ 



Si{m, r, r') 



m 



r+r' 



dyPi 



^2 _|_ y2, 

2rV 



2rr' J ' 
2£ + l'^< ' 



(30a) 

(30b) 

(30c) 
(30d) 



where r< = min[r, r'] and r> = max[r, r']. The kernel Si{m, r, r') is actually the 
Green function of the nonrelativistic kinetic energy operator and takes a simple 
form while the kernel Qe{m, r, r') can be calculated analytically for each value 
of i [37,38]. We have, here also, introduced in (30b) the parameter a which 
takes the value 1 respectively 2 for one respectively two (identical) particle 
problems. We can use an additional approximation and majorize Ki{y) by 
1/y. In this case, Ti{m,r,r') is replaced by A(_{m,r,r'), see (21b). Note that 
T^(0, r, r') — A^(0, r, r') since in this limit (m — > 0) mKi{my) — 1/y. 

As indicated in the previous Section 3.1, simplifications occur in the ultra- 
relativistic case (m = 0). With the help of the Holder inequahty we obtain 



Ne< B{n,e,p,p') 



Jo 



\V-{r)\ 




r dr\V-(r)\PP' 






Jo 



r 

Jo 



i/p 



drrP-^ \V- 



i/{pp') 

{n-3)/p 



(31a) 



with 



B{nJ,p,p') = [C{e,p/{p-l))r-'C{i,pp'/{p' - 1)), (31b) 

and where n > 2, p > 1, p' > 1 and C{i, q) defined by (22c). The upper limit 
on the total number of bound states that can be obtained from the upper 
limit on the number of £-wave bound states (31) is very similar to the upper 
limit (22a) except that the sum would stop aX ^ = L (the largest value of 
the angular momentum I for which bound states do exist) and that n is only 
restricted to be greater than two. 

Taking the limit n — > oo, we obtain the following necessary condition for the 
existence of at least one £-wave bound state: 



dr 



r 



C{l,p/{p-l))r\V-{r) 



IP 



> 1. 



(32) 
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Table 1 

Values of the function c{i), see (34), for several values of £. Note that c(0) = oo and 
c(oo) = \/27r. 



£ 


c(i) 


£ 


c(£) 


£ 




£ 




I 




s 






2 5S1 


40 


7 51 5 

. (J J. (J 


2 


2.795 


9 


2.550 


16 


2.529 


50 


2.513 


3 


2.678 


10 


2.545 


17 


2.528 


60 


2.512 


4 


2.625 


11 


2.541 


18 


2.526 


70 


2.511 


5 


2.596 


12 


2.538 


19 


2.525 


80 


2.511 


6 


2.578 


13 


2.535 


20 


2.524 


90 


2.510 


7 


2.566 


14 


2.533 


30 


2.518 


100 


2.510 



This simple relation yields a lower limit, g^^ on the (numerically) exact "criti- 
cal" value of gc,ex, for which a first bound state appears, gc,cx > fi'c o-nd yields 
an upper limit L"*" on L. We can obtain a simpler expression by considering 
now the limit p — > oo. We find that a necessary condition for the existence of 
at least one ^-wave bound state is also given by 



C{e, 1)M>1, (33a) 

where 

M = max[r\V-{r)\]. (33b) 

This last necessary condition have only a sens for £ > because C(0, 1) di- 
verges. Of course, this does not mean that there always exist bound states for 
i = 0. For this value of the angular momentum, one needs to use the relation 
(32) to draw conclusions about existence of bound states. The constant C{£, 1) 
can be rewritten as 



c<e 1) = J- + <0 ,34) 

Qv^r(« + 3/2)^27+1' ^ ' 



The interest of this rewritten is that the function c{£) varies very slowly with 
i and can then be easily and usefully tabulated, see Table 1. A possible simple 
majorization of c{£) is c{£) < v^27r + ai'^ with a = 0.7 and b = —1.18 (valid 
at least for 1 < i < 100). Note also that the necessary condition (33) is 
the ultrarelativistic counterpart of the well known nonrelativistic necessary 
condition {£ + l/2)~^max[r^ |2mV^~(r)|] > 1 (see for example [40]). 
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3.3 Upper limit on the total number of hound states for central potentials in 
the ultrarelativistic regime 



In this section, we use the upper hmit (31) obtained in the previous section 
to infer a limit on the total number of bound states applicable only to central 
potentials in the ultrarelativistic regime. We choose in (31) n = 2, p = 2 and 
p' — > cxD and with the help of (23a) (which is certainly correct up to £ = 100) 
we obtain 

with 

M^max[\V-{r)\], (35b) 



roo 



drr\V-{r)\. (35c) 

To some extend, the limit (35) is the ultrarelativistic counterpart of the non- 
relativistic Bargmann-Schwinger upper bound [19,21]. The upper limit on N 
is obtained with a summation over the right-hand side of (35a) from £ = 
to £ = L taking into account the multiplicity of each £-wave bound states. To 
this end, an upper limit on L is needed. The best limit, which behaves linearly 
with the strength of potential, is obtained, not from (35a), but instead from 
the simple relation (33). We have 

N<^{L+ + 1)MI. (36) 

The optimal value for L+ is obtained by solving (33) but a neater, if generally 
less stringent, upper limit L'^'^ (L+ < L"^"^) is given by 

L++ = l(VrT8^-l), (37a) 

with 

c(l) 

s = ^M. (37b) 

The expression (37) for the upper limit on L is obtained using (33), (34) and 
the second part of the inequalities \/x + 1 > T[x + 3/2]/r[a; + 1] > y/x for 
X > 0. c(l) could even be replaced by c(L++) in (37b). The transcendental 
equation obtained is solved very quickly after few iterations thanks to the slow 
variation of c{i) as a function of i. If we believe that the relation (23a) is always 
true for all values of the angular momentum, we can write the asymptotic 
expression of the upper limit (36) when the strength of the potential, g, goes 
to infinity: 

N{g ^ oo) < ^MMI, (38) 

where c{oo) were used to obtain the asymptotic expression of L"*", L+ pa A^/a, 
where the symbol ~ means asymptotic equality. This last expression is useful 
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to compare the upper limit (36) with the Daubechies upper hmit (24) since 
we clearly have 

/ drr'^\V-ir)\^ <MMX. (39) 

This last inequality means that the upper limit (24) would always be better 
than the limit (36) if its coefficient was equal to 1. Since the coefficient is 
greater than unity, A-kK = 1.294 > 1, there is some room for the limit (36) to 
be more stringent. A square well potential is an example: 

V(r) = -Vo e[{Ri - r)(r - R2)], (40) 

where R2 > Ri are two arbitrary positive radius and 9{x) is the step function 
{0{x) = 1 if X > and 6{x) = if x < 0). When the ratio of the radius, 
R1/R2, is in the interval [0.4859, 1), the upper limit (36) is better than the 
Daubechies upper limit. 



4 Tests 

In this section, we test the limits obtained in the previous sections with two 
simple central potentials written hereafter as V{r) — —gR~^v{r/R). This first 
one is an exponential potential 

v{x) — exp{—x). (41) 

The second one is a Poschl- Teller potential 

vix) = — ^ 
cosh X 

The upper limits (21) and (22) on the total number of bound states behave 
with the strength of the potential, g, at least as and should not be very 
stringent for strong potential possessing many bound states. Nevertheless, for 
weak potentials these upper bounds could be rather effective, especially to 
estabhsh necessary condition for the existence of at least one bound state. 
The same remark applies to the upper limits (30) and (31) since they behave 
at least as g^ instead of g. 

4-1 The S-wave case 

For simplicity, we restrict our attention to the n = 2 case to test the upper 
hmit (30) in the S-wave case. This hmit, for £ — 0, reads 

11 



(42) 



Table 2 

Comparison, for the exponential and the Poschl- Teller potentials, between the crit- 
ical values, of 5 yielded by the limits (43), (44), the critical values, gc,D, of g 
yielded by the Daubechies upper limit (24) and the exact critical value gc,ex, ob- 
tained by solving numerically the spinless Salpeter equation. 



/3 


Exponential 


Poschl-Teller 




9c,D 


9c,ex 


9c 


9c,B 


9c,cx 





4.443 


4.370 


5.574 


4.126 


3.886 


5.008 


1 


1.223 


0.6574 


1.361 


1.512 


0.8631 


1.742 


2 


0.6739 


0.3374 


0.7133 


0.8912 


0.4582 


0.9598 


3 


0.4604 


0.2261 


0.4804 


0.6233 


0.3092 


0.6549 


4 


0.3487 


0.1698 


0.3616 


0.4769 


0.2329 


0.4956 


5 


0.2803 


0.1360 


0.2898 


0.3854 


0.1867 


0.3981 



1 roo foo r- 

No<——g^ dxv{x) dyv{y) {Ko{(3\x-y\)-Ko{(3{x + y))) 
(ttq;)^ Jo Jo l 



+ t:(3{x + y — \x — y\ 



(43) 



with /? = mR. g and (3 are two dimensionless parameters. While the depen- 
dency of (43) on g is very simple, the behavior of this limit with a variation 
of /? is more involved (except for large /3, see below). Simplifications obviously 
occur in the case of a vanishing mass and the limit reads 



Nn < 



[nay 



dxv(x) 



dyv{y) log^ 



X + y 



x-y 



(44) 



We test the above limits by computing the necessary condition on g, g > gc, 
for the existence of at least one bound state. In other words, we calculate 
a lower limit, g^, on the (numerically) exact "critical" value of g, gc^cx, for 
which a first bound state appears, g^cex > Oc- The numerical calculations are 
performed with the help of the very accurate Lagrange Mesh method [41]. The 
comparison between the (numerically) exact results, ^'cex, and the lower hmits, 
gc, is displayed in Table 2 for a two-particles problem (a — 2). We choose a 
two identical particles problem only for numerical convenience (to not modify 
numerical codes). The results obtained with the upper limits (43) and (44) 
are rather satisfactory especially when /3 is large and are always better than 
the results obtained with the Daubechies upper hmit (24). 

In the limit of large f3, only the kernel iSo(m, r, r') will obviously contribute 
significantly. In this asymptotic regime, the lower limit g^ takes a very simple 
form, and reads for the exponential (E) and the Poschl-Teller (PT) potentials 
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E \/2 1.4142 , PT 1.9663 
9.^«^ = — and sr«— (45) 

where the symbol ~ means asymptotic equality. But the limit of large /3 is 
equivalent to the limit of large m, which leads to nonrelativistic regime. We can 
thus compare the lower hmits (45) with those obtained from nonrelativistic 
upper bounds on the number of bound states, and also with the nonrelativistic 
exact results. The most stringent nonrelativistic lower bound QcNR is obtained 
with the upper limit found in [27] (see also [32]) and reads 

p 1.4383 , PT 1.9910 , 
^c,NR = and ^,NR = (^^) 

The exact result, for the exponential potential, is obtained in terms of the first 
zero, Zo = 2.4048, of the Bessel function Jo (a;) (see for example [42, p. 196]) 
while the exact result, for the Poschl-Teller potential, is extracted from the 
exactly known spectra (see for example [42, p. 94]). These exact results read 

E _ 4 _ 1-4458 1 PT _ 2 

5'c,NR,ex ~ 4^ ~ p ^ 5'c,NR,ex " ^- (4') 

At first sight, it is surprising that the upper limit on the number of bound 
states derived for spinless Salpeter equation (which reduces to the Schrodinger 
equation in the limit of large m) leads asymptotically to similar restrictions 
on g, for the existence of at least one bound state, than previously known 
nonrelativistic upper limits. It should be noted that the limits obtained in this 
article has no counter part in the nonrelativistic context since we have used 
iterated kernels to obtain these upper limits, and in particular the limit (30). 
Indeed, it is not necessary to use iterated kernels to obtain an upper limit on 
the number of £-wave bound states in the context of nonrelativistic quantum 
mechanics (the Bargmann-Schwinger upper limit [19,21]) since the kernel is 
not singular. Moreover, the use of an unnecessary iterated kernel would lead 
to an even more incorrect dependency on g (the Bargmann-Schwinger upper 
limit docs not display the correct dependency on g, see for example [22]), which 
would lead to a poorer upper limit for a strong potential possessing many 
bound states. Nevertheless, it seems that such a poor upper limit could yield 
very strong conditions on g for the absence of bound states. This remark and 
the potential benefit which could be obtained from it, namely obtaining strong 
conditions for the absence of bound states in the context of nonrelativistic 
quantum mechanics, will be studied in detail elsewhere. 

These comparisons with nonrelativistic results and those reported in the Ta- 
ble 2 show that the upper limit (30) is rather effective to establish necessary 
condition for the existence of at least one bound state for the spinless Salpeter 
equation, especially when (3 becomes large. This also indicates that this upper 
limit should yield cogent information for weak potentials possessing few bound 
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Table 3 

Comparison, for the exponential and the Poschl- Teller potentials, between the crit- 
ical values, ^c, of g yielded by the limits (32) {gc,i) and (33) (^cii) and the exact 
critical value 5c,cx, obtained by solving numerically the spinless Salpeter equation. 



I 


Exponential 


Poschl-Teller 


9c,l 


50,11 


5c, cx 


del 


50,11 


5c, cx 


1 


8.524 


6.922 


10.98 


7.437 


5.687 


9.545 


2 


13.67 


12.81 


16.39 


11.59 


10.53 


14.04 


3 


19.03 


18.46 


21.81 


15.91 


15.17 


18.52 


4 


24.44 


24.02 


27.24 


20.30 


19.73 


22.99 


5 


29.88 


29.53 


32.67 


24.73 


24.27 


27.46 



states. The results presented in Table 2, for n = 2, should be improved as n 
goes to infinity and, we beheve, should converge to the exact result if there 
were not the majorization (19). 

We can also test the upper limit (31) and the necessary condition (32) (both 
valid only when m = 0) by the computation of Qc- The series of Qc obtained 
with (31) for each vahic of n arc monotonic series which converges to (^c = 4 
for the exponential potential and converges to = 3.685 for the Poschl-Teller 
potential, which are, obviously, the values obtained with the expression (32). 

4-2 The £-wave case 

In this section we only test the relations (32) and (33) by computing the lower 
limit Qc for £ > with the potentials (41) and (42). The comparison of these 
lower limit with the (numerically) exact critical coupling constant is given in 
the Table 3. 

The lower limits yielded by the relation (32) arc always better than those 
obtained with the relation (33) but the differences become smaller as £ grows. 
These lower limits are quite satisfactory compared to the exact results and 
the relative differences between these quantities decrease from 22% to 10% for 
£ increasing from 1 to 5 for both potentials. 

4-3 Tests of an upper limit on L 

To conclude this section devoted to test the limits obtained in this article, we 
compare the upper limit, L"*", obtained with (33), with the exact largest value, 
L, of the angular momentum for which bound states do exist. The results 
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Table 4 

Comparison, for the exponential and the Poschl- Teller potentials, between the upper 
limit, L+, obtained with (33), with the exact largest value, L, (obtained by solving 
numerically the spinless Salpeter equation) of the angular momentum for which 
bound states do exist. 





Exponential 


Poschl- Teller 


9 




L 


L+ 


L 


10 


1 





1 


1 


20 


3 


2 


4 


3 


30 


5 


4 


6 


5 


40 


6 


6 


8 


7 


50 


8 


8 


10 


10 


100 


17 


17 


21 


21 


150 


27 


26 


33 


32 


200 


36 


35 


44 


43 



obtained for an exponential and a Poschl- Teller potential are reported in the 
Table 4. The bounds on L obtained with (33) are very stringent for both 
potentials. These excellent results are not so surprising since the same strong 
limitations on L are obtained with the nonrelativistic counterpart of (33). 



5 Upper limit on the number of bound states lying below a given 
energy 



To conclude this letter, it is worth mentioning that the method used to obtain 
the upper limits presented in this work can also be used to derive upper limits 
on the number of states of a given potential that lie below a given energy. 
Indeed, instead of using — |V^~(?")| as a comparison potential, which admit at 
least as much bound states as V(r), we can use a potential V{r) that equals 
V{r) wherever V{r) < —k^ {k^ is an arbitrary real energy) and equals —k^ in 
the regions for which V{r) > —h?. But this latter potential, V{r), is equivalent 
to the potential Vi^{r) that equals V{r) + wherever this quantity is negative 
and equals zero otherwise, since V^ir) = V{r) + k^. Thus, an upper bound to 
the total number of £-wave bound states associated with V^{r) serves to limit 
the number of £-wave states that lie at or below the energy E — —k^ for the 
potential V{r). For example, the limit (30) for n = 2 reads 



rco poo 

Ni{E<-^)< drV,{r) dr' V^{r')Tf{ 
Jo Jo 



m, r, r 



(48) 

V+K^<0 
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with T£{m,r,r') defined by (30b). This last upper hmit can be tested with 
the harmonic oscillator potential, V{r) = k'^'r^, with m = (with a slight 
adaptation since V{r) is now positive). Indeed, in this case the spectrum is 
exactly known ior i — [43] : 

En^kXr,, (49) 

where — A„ are the zeros of the Airy function. Thus, below these eigenenergies, 
there are exactly n — 1 bound states. For n — 1,2 and 3, the upper limit gives 
1, 8, 21. As expected, the limit gives cogent results when few bound states 
exist and becomes less stringent as the number of bound states grows because 
of the incorrect behavior with the growth of the coupling constant. 
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